Thermally assisted adiabatic quantum computation 
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We study the effect of a thermal environment on adiabatic quantum computation using the Bloch- 
Redfield formahsm. We show that in certain cases the environment can enhance the performance in 
two different ways: (i) by introducing a time scale for thermal mixing near the anticrossing that is 
smaller than the adiabatic time scale, and (ii) by relaxation after the anticrossing. The former can 
enhance the scaling of computation when the environment is superohmic, while the latter can only 
provide a prefactor enhancement. We apply our method to the case of adiabatic Grover search and 
show that performance better than classical is possible with a superohmic environment, with no a 
priori knowledge of the energy spectrum. 



I 

O 

o 



> 

m 
m 

o 

o 



I 

O 

o 



Quantum computation (QC) aims to harness the phys- 
ical resources made available by quantum mechanics to 
gain an advantage over classical computation. A major 
obstacle to construction of a large scale quantum com- 
puter is loss of coherence resulting from uncontrolled cou- 
pling to the environment. In principle, environmental ef- 
fects may be circumvented by the use of quantum error 
correction [J, 0, [3, 0] ■ In practice, however, such schemes 
require significant overhead. It is therefore likely that 
many noisy qubits will be available before many error- 
corrected logical qubits are. 

This observation motivates the search for models of 
QC with intrinsic resistance to noise. One such example 
is adiabatic quantum computation (AQC) [1, @, Here 
we investigate a regime in which weak coupling to an 
environment can improve the performance of AQC. 

In AQC, information is stored in the ground state of 
a quantum system and manipulated by control of the 
system Hamiltonian. An AQC is operated by deforming 
an initial Hamiltonian Hi into a final Hamiltonian Hf 
through intermediates Hs = [1 — A(i)]7?i -I- X{t)Hf, with 
A(t) changing from to 1 between the initial {ti=0) and 
final (tf) times. If the evolution satisfies the adiabatic 
condition {h=kB=l throughout): \{l\dH / dt\{))\ < .g^(A), 
where 5(A) is the energy gap between the ground (|0)) 
and first excited (|1)) states, then the system will be in 
the ground state of Hf at tf with probability close to 



one, and the solution may then be read out [2l|. In a 
global adiabatic scheme, A — t/tf and the adiabatic con- 
dition must be satisfied for the smallest gap gm- If 5(A) 
is known, one can choose dX/dt oc g^{X) to enhance the 
performance using a local adiabatic scheme Here, we 
assume no a priori knowledge of the energy spectrum, 
and use A = t/tf throughout. The amount of time re- 
quired to successfully run a computation is determined 
by the minimum gap between the first two energy levels, 
gm, along the path connecting Hi and Hj. In order for 
the evolution to remain adiabatic throughout, the total 
time required is tf (x !/(?„. 

Here we analyze the behavior of an AQC in the pres- 
ence of a thermal environment with temperature T ^ 
g-m- We restrict our analysis to problems in which the 



performance is limited by a single minimum gap of the 
type of an energy level avoided crossing. This corre- 
sponds to a first order quantum phase transition, which 
is believed to be hardest for AQC [lo| . 

In general, if there are I energy levels within the range 
T from the ground state, then thermalization can sup- 
press the ground state probability by at most a factor of 
l^^. For a Gaussian distribution of the levels, I is polyno- 
mial in the number of qubits n, if T is much smaller than 
the total spectral width. In this case, one may compen- 
sate for thermalization by repetition with a polynomial 
overhead. Moreover, the transition times are expected to 
be very long, probably longer than the computation time, 
otherwise classical annealing would yield the solution ef- 
ficiently. This is different from an anticrossing for which, 
as we shall see, the transition rate is sharply peaked at 
the anticrossing point. We therefore only focus on the 
anticrossing and use 2-level approximation. 

Let us assume that the minimum gap occurs at A=A„i. 
We adopt a new coordinate, e=2i?(A— A„i), where E is an 
energy scale characterizing the anticrossing. Close to the 
anticrossing, the system Hamiltonian within the 2-level 
approximation is well described by 
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and the gap between the first two states is well approxi- 
mated hy g — \/ + g^ ■ Here Tx,z are the Pauli matrices 
in the 2-level subspace. Due to the Landau-Zener transi- 
tion [III, [l2| , the probability of being in the excited state 
at t = tf is given by 
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where ta = 4:E/TTg^ is the adiabatic time scale (see Table 
Ufor definition of all time scales). 

We incorporate the environment by assuming that 
qubits are coupled to bosonic heat baths that are in equi- 
librium with gm T <^ E. The total Hamiltonian is 
H = Hs + Hb + Hint , where Hb and iJint are bath and 
interaction Hamiltonians respectively. We also assume 
that in the 2-level subspace the interaction Hamiltonian 
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TABLE I: Characteristic time scales for t 



f 



ta time scale imposed by non-adiabatic transitions 

td time scale imposed by thermal mixing 

tr time scale imposed by relaxation after the anticrossing 



has the form 



(3) 



where Q is an operator representing the collective effect 
of all baths on the 2-state problem. Equations ([T]) and 
dS]) capture the physics of a wide range of problems that 
have one sharp anticrossing. 

For slow evolutions of the Hamiltonian considered here, 
as long as the correlation time of the environment is 
shorter than decay times of the s yste m, one can safely 
assume Markovian approximation |22l | . Writing the den- 
sity matrix as p = (l + p-'r)/2, the 2-state Bloch-Redfield 
equations are [l4| : 
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where 7 = {9ml9?[S{9) + S{-g)l 7^ = ^{^l9?S{Q), 
7 = 7 + 7^, and peq ^ [S{9) - S{-g)]l[S{g) + S{-g)]. 
Here, the bath's spectral density is defined as S{uj) = 
dt e"^^{Q{t)Q{0)), where (...) denotes averaging over 
environmental degrees of freedom. The prefactor {gm/ gY 
makes 7 sharply peaked at e = as expected. 

For a bosonic environment fl3|, S{uj) = J{uj)/{1 — 
e"'^/"^), where J{uj) = rjujluj/iVcl'^ e^^^^^, with lOc being 
a cutoff frequency which is assumed to be larger than all 
other relevant energy scales in the system. Therefore 

7 = i9m/gfJig) coth(5/2r), p^q = tanh(g/2r). (5) 

Here, we focus only on ohmic (s = 0) and superohmic 
(s > 0) cases for which the correlation time of the bath 
(~ 1/ujc) is short compared to the relevant time scales. 
A subohmic (s < 0) environment has a large correlation 
time, hence the Markovian approximation and therefore 
Bloch-Redfield equation do not hold 

We are interested in problems with small gap, 
gm'^T<^E. We divide the evolution into three regions, 
as shown in Fig. [TJ In region I, the gap is larger than 
T and thermal transitions are suppressed. In region II, 
both thermal and non-adiabatic transitions between the 
two states are possible. In region HI, the system again 
has a gap larger than T, but now the system can relax 
from the excited state to the ground state. Such relax- 
ation can only increase the probability of success. 

Let us start by finding the excitation probability im- 
mediately after region II. Assuming T^g, which holds 




FIG. 1: Ground and first excited states of a system with one 
anticrossing at A = 0.5. The 3 regions are identified in such a 
way that g > T for regions I and III, and g < T for region II. 



for most of the region, we have p^q — 0, and 
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where "fo^2r]T. We perform the calculation in the regime 
gm,"^!- The presence of the damping terms in the first 
two equations in ([4]) will make px and py decay in a time 
scale 1/7) much shorter than the relevant time scale 
for Pz ^/gm)- Thus, to find the slow evolution of pz, 
one can use the stationary values for p^ and py, obtained 
from px^y = 0. Solving the first two equations in ^ for 
Py and substituting into the third equation, we get 
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Here, F is the rate of transition between the two branches 
of energy that meet at the anticrossing. Integrating ([7]), 
we find 

,„Ml = -r.<.r(„.-/""irw, (8) 

Pz{0) Jo J-T e 



which leads to 



p(e « T) = p(0)e~*//*^ 



(9) 



where is a characteristic time scale that describes the 
thermal mixing near the anticrossing, and is given by 
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We have taken the integration limits to infinity assuming 
that F(e) is sharply peaked at e = 0, which is the case 
for s < 1. Using the initial condition pz{0) = 1, the 
excitation probability after region II is approximately 



Pi(e«r) = i(l + e-*^/*-). 



(11) 



For a linear time evolution (i.e., A = t/tf), one has 
e = 2E/tf. If the environment is ohmic, then s = and 
7 = 7o is a constant. Therefore 
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in agreement with Kayanuma [l3| and Ao and Rammer 
[ist . For fast evolutions (short tf), (|lip behaves the same 
way as a closed system ([2]) , while in the slow regime (long 
tf), PT|) gives a ground state probability of ~ 1/2, cor- 
responding to the complete mixture of the two states. 

For a superohmic bath (s > 0), = 0, hence J — J — 
lf)(.9m/gY{g/^cy- In such a case. 
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where A = jog'^/uj^. The important contribution 
to the integral comes from regions with g « |e| ^ 

Ai/(3--) » 

gm, where the inequality follows from gm < 
lo{gm/i^cY = 7(0), which was our initial assumption. 
This condition, however, can be satisfied in the limit of 
9m 0, only if s < 1. Replacing e — A^/^'^^^^x, we find 
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Since as is independent of gm, we have td cx gm^'^'^ 
which scales better than ta ~ g^ ■ It is easy to check 
that the integral is convergent for s < 1. For s > 1, 
the condition g„i <C 7 cannot be satisfied in the limit of 
g-m, — ^ 0, invalidating our approach. 

We now study the effect of relaxation after the anti- 
crossing (region III). From ([7]), we see that r(e^7) « 7. 
The probability of ending up in the excited state be- 
comes Pi/(e w £;) w Pi(e w T)e^^T Using ^ 
and assuming coth((7/2T) « 1, which holds for most of 
the region, and e 3> .g™, we find 
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Here we have defined a third time scale tr that character- 
izes such a relaxation process. One can write tr=ta/ Ks, 
where 
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\ti{E/T) s = 



(16) 



Notice that tr slowly decreases with T. 

The probability of success, i.e., the final ground state 
probability, is therefore given by 
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It reaches ~ 1/2 in a time tf ^ td, but approaches 1 
in a time tf ^ tr. If id < 7 it is advantageous to run 
the system faster but repeat the process. The relevant 
time scale for computation will then be ~ 2td, which for 
an ohmic environment is ^ ta (x l/gm, the time scale 
for a closed AQC. On the other hand, for superohmic 
cases with 0<s<l, one has td cx gm^^^^ which shows 
an improved performance compared to the closed AQC, 



as gm 0. The performance becomes better as s gets 
closer to zero, until s = (i.e., ohmic) at which point 
the low frequency part of the noise spectrum becomes 
nonzero and the performance goes back to 1/ gm- 

This sudden change at s = is related to the sharp 
jump in 5(0;) oc at o; = 0, from a nonzero value 
at s = to zero at s > 0. However, the S'(O) that 
appears in the definition of 7^ is not exactly zero fre- 
quency, but really the low frequency component of the 
noise, i.e., S'(~l/i/). As s becomes smaller, the low fre- 
quency component gets larger and eventually dominates 
the 7 in ([8]), resulting in a smooth transition to the ohmic 
behavior. Without the 15(0) term, an ohmic environment 
would yield a ~ gm^^ behavior. Here, a competition 
between pure relaxation, which tends to enhance the per- 
formance, and pure dephasing (due to the low frequency 
noise) which works against it, is noticeable. Taking both 
processes into account, in the case of ohmic environment, 
the performance of the system will be the same as that 
for a fully coherent AQC. 

For systems with td > tr, the computation time scale 
will be determined by tr and the ground state probability, 
for small tf, will basically have the form Pofit) « 1 — 
g-t//tr^ If tr < ta, then we will again have a better 
performance compared to a closed AQC. However, as we 
saw before, tr has the same dependence as ta- Thus, 
any speedup over AQC by this process can only be via 
a prefactor Kg (if it is larger than 1). The enhancement 
reported in Ref. ^ falls in this category since the number 
of qubits considered was not large enough to obtain small 
gm and therefore thermal mixture at the anticrossing. 

We should emphasize that equation (fT6|) is calculated 
assuming that the 2-state approximation holds for the 
entire range. While this can be the case for some Hamil- 
tonians, such as adiabatic Grover search , it is not true 
in general. In fact, it is very difficult to calculate tr for 
a general problem. However, one would not expect this 
type of relaxation, which is equivalent to classical anneal- 
ing, give any scaling benefit over classical computation. 

We now apply our approach to the adiabatic implemen- 
tation of Grover's search algorithm 0, In this case, 
the explicit dependence of g„i on the problem size may 
be obtained, and hence all quantities may be calculated 
in terms of the size of the unstructured search problem 
N = 2". Following Roland and Cerf [|, we use the 
Hamiltonian Hs = £; [1 - (1 - A) \+) {+\ - A |m) (m|] , 
where |m) is the marked state to be found and 
\+)={N)'^/^J2i\l)- Defining e=£;(2A - 1), the gap 
is g{e)=^E^/N + (1 - 1/A^)e2. The minimum gap. 



gm—E/vN, hcs at e = 0. The third energy level, E2—E, 
has (A^— 2)-fold degeneracy. Aglobal adiabatic algorithm 
{X=t/tf) results in tf^N/E [7|. Using a local adiabatic 
algorithm 0, one can achieve tf^VN /E. Because of the 
large degeneracy of E2, the 2-level approximation will 
only be valid in the temperature regime T<^E/ log A^. 
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FIG. 2; Probability of success for AGS with 12 (blue), 16 
(red), and 20 (green) qubits, and a superohmic environment 
with s=0.5 at r/i?=0.1. The solid and dotted curves are for 
longitudinal (77^=0, f;^=0.1) and transverse (jf=Q.\, 'ff=^) 
couplings to the environment, respectively. The solid black 
line represents a closed system (ff=ff=Q'), which with the 
normalized x-axis is independent of n. 



We consider the implementation of the unstructured 
search problem on n qubits, and hence may describe our 
noise model in terms of operators acting on these qubits. 
We show that the type of 2-level noise model (3) indeed 
arises for a general coupling of qubits to the environment: 
-ffint ^ ~Y^=A^i® ^Zi^al], where trf are the Pauli 
matrices for the i-th qubit, and X^, Zi are its correspond- 
ing heat bath operators. In the large N limit, the effec- 
tive 2-level system and interaction Hamiltonians become 
([H) and ([3]), respectively, where Q — \ — Zi). As- 

suming uncorrelated heat baths, ([T7]) also holds for this 
problem with 77 = \n(ff ^ ff \ where ff'^ are average 
friction coefficients for the Xi and Zi operators. 

For large n, the scaling of td with N is given by id iV 
for ohmic, and ~ 7V2/(3-s) ^^j. g^perohmic environ- 
ment (with linear interpolation). It is clear that for su- 
perohmic environment with s < 1, the scaling is better 
than that for classical computation. 

We have also performed numerical simulations of adia- 
batic Grover search, solving the Bloch-Redfield equations 
without the 2-level or large TV approximation for 12, 16 
and 20 qubits. Figure [2] plots Pq/ as a function of i//ia 
for a case with superohmic environment with s = 0.5. As 
is clear from the figure, the curves increase faster com- 
pared to a closed system for larger n (smaller Qm). This 
agrees with the scaling advantage of the noisy system 
compared to the closed system according to our analyti- 
cal prediction. 

To summarize, using the Bloch-Redfield formalism we 
have identified 3 time scales for the evolution of AQC 
and determined their scalings with gm- We have shown 
that relaxation after the anticrossing can only provide a 
prefactor enhancement for computation time. Thermal 
mixing at the anticrossing, on the other hand, can en- 
hance the scaling of the computation if the environment 
is superohmic with 0<s<l, while the same environment 



will be destructive for gate model QC. This underlines 
the important difference between the two models in re- 
sponse to the environment. Finally, we should mention 
that a presence of low frequency noise, as in spin envi- 
ronment [20| . will remove the above enhancement. 

Discussions with A.J. Berkley, J.B. Biamonte, A. Blais, 
E. Farhi, E. Ladizinsky, A.J. Leggett, D. Lidar, A. 
Maassen van den Brink, D.A. Meyer, G. Rose, A.Yu. 
Smirnov, P.C.E. Stamp, and M. Wubs are gratefully ac- 
knowledged. MHSA would also like to thank D.V. Averin 
for pointing out the importance of low frequency noise 
and the breakdown of the secular approximation in the 
Bloch-Redfield formalism. 
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